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Abstract: Results of inclusive measurements of inelastic electron and positron scattering 
from unpolarized protons and deuterons at the HERMES experiment are presented. The 
structure functions -Fg ^'^d F2 are determined using a parameterization of existing data for 
the longitudinal-to-transverse virtual-photon absorption cross-section ratio. The HERMES 
results provide data in the ranges 0.006 < X < 0.9 and 0.1 GeV^ < < 20 GeV^, covermg 
the transition region between the perturbative and the non-perturbative regimes of QCD 
in a so-far largely unexplored kinematic region. They are in agreement with existing world 
data in the region of overlap. The measured cross sections are used, in combination with 
data from other experiments, to perform fits to the photon- nucleon cross section using the 
functional form of the ALLM model. The deuteron-to-proton cross-section ratio is also 
determined. 
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1. Introduction 



Over the past decades, lepton-nucleon scattering has played a major role in the development 
of our present understanding of nucleon structure. For a review on the subject see for 
example In lowest order perturbation theory, scattering of charged leptons I off nucleons 
N proceeds via the exchange of a neutral boson (7*, Z^). At the HERMES lepton-nucleon 
centre-of-mass energy of ^/s = 7.2 GeV, contributions from Z*^-exchange to the cross section 
can be neglected. Therefore, only the electromagnetic interaction in the approximation 
of one-photon exchange is considered here. In this approximation, the differential cross 
section of unpolarized inclusive charged-lepton-nucleon scattering, I + N ^ I' + X (where 
X denotes the undetected final state), is parameterized by two structure functions Fi{x, Q^) 
and F2{x,Q'^). Here x = Q'^/2Mv is the Bjorken variable, with —Q^ being the square of 
the four-momentum transferred by the virtual photon and v its energy in the target rest 
frame. The variable a; is a measure for the inelasticity of the process with < x < 1, and 
X = 1 for elastic scattering. 

In the deep-inelastic scattering (DIS) regime, \/Q^ and are much larger than the 
typical hadronic scale, usually set to be the mass M of the nucleon, and the invariant mass 
W of the photon-nucleon system is much larger than the masses of nucleon resonances. In 
the Quark-Parton Model (QPM), the DIS process is viewed as the incoherent superposition 
of elastic lepton scattering from quasi-free point-like quarks of any flavor q. The variable 
X can then be interpreted as the fraction of the longitudinal nucleon momentum carried 
by the struck quark in a frame where the nucleon moves with infinite momentum in the 
direction opposite to that of the virtual photon. In this picture, quark distribution functions 
fq{x,Q'^) describe the number density of quarks of flavor g in a fast-moving nucleon at a 
given value of (x, Q^) and experimental values of F2{x,Q'^) have been used to constrain 
these. At low values of Q^, where this picture of incoherent quasi-free scattering does not 
apply, phenomenological models have been developed (see e.g. Refs. and |Q) to describe 
the measured structure functions. 



A wealth of unpolarized inclusive charged-lepton DIS data is available from the collider 
experiments HI g |, g 0, i and ZEUS §, ^, 0, |l3 at HERA with lepton- 

nucleon centre-of-mass energies ^/s up to 320 GeV, the muon experiments BCDMS [|l5|, 

31 GeV), experiments with electrons at 
I, m) (Vs < 3.25 GeV). The HERMES 



EMC ig, NMC and E665 |18| (^ ^ 12 
SLAG [||] {^/s<7 GeV) and at JLAB |21 



experiment [24| at HERA collected a large data set for positron and electron scattering 
on a variety of nuclear targets, including the proton and deuteron data presented here. In 
particular, the HERMES data cover the transition region between the perturbative and 
non-perturbative regimes of QCD in a kinematic region so far largely unexplored by other 
experiments. In this work, these data are presented together with fits to the world data for 
the photon-nucleon cross section using the Regge-motivated approach of the ALLM |^ ^ 
model. The paper is organized as follows: the formalism leading to the extraction of 
the structure function F2 is briefly reviewed in Sect. ^; Sect. ^ deals with the HERMES 
experimental arrangement and the data analysis is described in Sect. 0. The systematic 
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uncertainties in the resulting cross sections and cross-section ratios are discussed in Sect. |5|. 
Section ^ offers a discussion of the results and Sect. |^ provides a summary. 



2. Formalism 



mi 



M 



k = {E,k), k' = {E',k') 



P ^= (M, 0) 



q = k-k' 



lab 



y 



M 

Q' 

2P-q 2Mv 
P ■ q Uh f 



P-k E 

W"^ = {P + qf = M2 + 2Mv - Q2 



Lepton mass (taken to be negligible) 
Mass of target nucleon 

4-momenta of the initial and final state leptons 
Polar and azimuthal angle of the scattered lepton 
4-momentum of the initial target nucleon 
4-momentum of the virtual photon 

Negative squared 4-momentum transfer 

Energy of the virtual photon in the target rest 
frame 

Bjorken scaling variable 

Fractional energy of the virtual photon 

Squared invariant mass of the photon-nucleon 
system 



Table 1: Kinematic variables used in the description of lepton scattering. 

In the approximation of one-photon exchange, the inclusive differential cross section 
for scattering unpolarized charged leptons on unpolarized nucleons can be conveniently 
parameterized in terms of the structure functions Fi{x^Q'^) and F2{x,Q'^): 



dx dQ"^ 



(2.1) 



where aem is the fine-structure constant and all other variables are described in Tab. |l|. 
The quantities x and are fully determined by the kinematic conditions of the incident 
and scattered leptons and the target nucleon. This cross section can also be written in 
terms of longitudinal (L) and transverse (T) virtual-photon contributions 



T[<jT{x,Q^) + eaL{x,Q^)] , 



(2.2) 



dx dQ"^ 

where ai and (Tj- are the absorption cross sections for longitudinal and transverse virtual 
photons, r is the flux of transverse virtual photons and the virtual-photon polarization 
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parameter e is the ratio of virtual-photon fluxes for longitudinal and transverse polariza- 
tions [^]. The structure functions Fi{x,Q'^) and F2{x,Q^) can then be expressed in terms 
of the two virtual-photon absorption cross sections (Tl{x^Q'^) and aT{x,Q'^): 

Fi{x,Q^) = -^ MK-aT{x,Q^), (2.3) 

F2{x,Q^) = -^ '^^L_.[a^^^^Q^)+aT{x,Q^)] , (2.4) 

4^ Olem ^ ^ V 



4M2 



x^ 



where K = i/(l — x) in the Hand convention |2?, The longitudinal-to-transverse 



photon- absorption cross-section ratio R = a^/ ax can be expressed in terms of Fi and F2: 

A determination of the structure functions i^i(x, Q^) and F2{x, Q'^) requires in principle 
cross-section measurements made at the same x and but at two or more different values 
of y (see Eq. ( |2.lD ), i.e., with different beam energies. The HERMES data used for this 
analysis were taken at a single beam energy. In such a situation, it is common practice to 



re-parameterize the cross section, Eq. (|2.l|) , as a function of F2 and R using Eq. (|2.5D : 
(fa A7ral^F2{x,Q 



dx dQ^ 



Q2 ^ y' + QVE'^ 
1-y- + 



2.6) 



4£;2 2[l + i?(x,Q2) 

The structure function F2{x,Q'^) can then be extracted from a single cross-section mea- 
surement at a given (x,(5^), by using a parameterization for R{x,Q'^) obtained from the 
available world data. This approach has been used in the analysis presented in this paper. 

3. The Experiment 

The HERA facility at DESY comprised a proton and a lepton storage ring. HERMES was 
a fixed-target experiment using only the lepton beam, which consisted of either electrons 
or positrons at an energy of 27.6 GeV, while the proton beam passed through the non- 
instrumented horizontal mid-plane of the HERMES spectrometer. An open-ended storage 
cell that could be fed with either polarized or unpolarized gas was installed internally to 
the lepton ring. 

The HERMES spectrometer, which consisted of two identical halves above and below 
the electron beam, was a forward spectrometer [^] with multiple tracking stages before and 
after a 1.5 Tm dipole magnet. It had a geometrical acceptance of ±170 mrad horizontally 
and ±(40 — 140) mrad vertically for particles originating from the center of the target 
cell, resulting in polar scattering angles 9 ranging from about 40 to 220 mrad. Particle 
identification (PID) capabilities were provided by combining the responses of a lead-glass 
calorimeter, a pre-shower hodoscope (H2), a transition-radiation detector (TRD), and a 
threshold Cerenkov detector that was upgraded to a dual-radiator ring-imaging Cerenkov 



detector (RICH) ||2^, 3C] in the year 1998. The lead-glass calorimeter and the pre-shower 



hodoscope were included in the trigger together with two other hodoscopes (HO and HI). 
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In this experiment target gases of hydrogen and deuterium were used. Part of the data 
were taken with polarized hydrogen and deuterium, with the target spin being reversed in 1- 
3 min time intervals so that the target was effectively unpolarized. In the case of hydrogen, 
using areal densities of the order of 10^^ nucleons cm~-^ and lepton currents of typically 
about 30 mA, luminosities of the order of 2 • 10'^^ cm^^s^^ were achieved for the polarized 
running, and about 10 times higher values for unpolarized running. The luminosity was 
measured by scattering the lepton beam off the atomic electrons of the target gas, i.e., 
M0ller scattering e~e~ — )■ e~e~ for an electron beam and Bhabha scattering e~^e~ — >• e'^e~ 
together with the annihilation process e~^e~ — )• 77 for a positron beam. The cross sections 
for these processes are precisely known in Quantum Electrodynamics, including radiative 
corrections. The scattered particles were detected in coincidence by two identical small 



calorimeters |31] located symmetrically with respect to the beam pipe. The coincidence 
rate of the pairs of leptons (and photons) provided a relative monitor of the luminosity. 
An absolute calibration of the luminosity measurement was provided by correlating the 
coincidence rate with the yields of the M0ller, Bhabha and annihilation processes. 

4. Data Analysis 

An event is accepted if it contains a track identified as a lepton by the PID system (see 
Sect. [4.3| ), and satisfies the selection criteria described in Sect. The number of mea- 
sured events Nmeas in each (x,Q^) bin is corrected by subtracting the charge-symmetric 
background from secondary processes Ncs and by dividing the resulting number by the 
corresponding trigger and lepton-identification efficiencies Strigger and Siep, while taking 
into account the hadron contamination Chad- 

^events {-^meas ^cs) ' ~p, ' 'p, • (^'i^) 

^trigger ^lep 



These corrections are described in Sects. 4.2 to 4.4 



The experimental cross section is then obtained as the ratio of the number of events 
^events in cach (x,(5^) bin of widths Ax and AQ^, and the integrated luminosity L (see 



Sect. 14.51) : 

An unfolding procedure for disentangling instrumental and radiative effects from the 
measured cross section is then applied in order to obtain the Born cross sections <7^fj,„ = 
(jP4 ^ggg Sect. |4.6| ). The structure functions F2 and F2 are finally derived from the Born 
cross sections through Eq. ( |2.6| ) using the parameterization R = -Riggg [32|. Two more 



corrections related to detector geometry and alignment are discussed in Sects. 4.7 and 4 
4.1 Event selection 

The kinematic range of the events selected for this analysis is shown in Fig. |l|, together 
with the requirements imposed on the kinematic and geometrical variables. 
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-2 -1 

10 10 X ^ 



Figure 1: Binning in (x, Q^) used in the analysis and kinematic acceptance of events. The kine- 
matic region covered is hmited by the geometrical acceptance in 9 and constraints on y and W'^. 
The symbols mark the locations of the average values of {x, Q^) of each bin. The symbols A to F 
denote bins with increasing Q2 at given X. 



Year 


events, 


in million 




proton 


deuteron 


1996 


2.3 


2.8 


1997 


4.6 


3.1 


2000 


9.5 


12.6 


total 


16.4 


18.5 



The tracks are required to be fully contained within 
the fiducial geometric acceptance of the HERMES spec- 
trometer. The constraint W"^ > 5 GeV^ excludes the 
region of nucleon resonances and acts as a selection 
for y > 0.1. The constraint y < 0.85 discards the 
low-momentum region, where radiative effects increase 
and where the trigger efficiency has not yet reached its 
plateau. The requirements imposed on y select the 
momentum range 4.1GeV< p < 24.8 GeV for the de- 
tected particle. The resulting {x,Q^) region, 0.006 < 
X < 0.9 and 0.1 GeV^ < < 20 GeV^, is subdivided 
into 19 bins in x and each x bin into up to six bins in 

Table |2| shows the numbers of events for each year 
of measurements used in this analysis, before the application of any of the corrections 
discussed in the following sections. 



Table 2: Number of raw events 
Nmeas uscd in this analysis, sepa- 
rated into the years of data taking. 
The numbers correspond to the total 
luminosities of about 450 pb~^ on the 
proton and about 460 pb~^ on the 
deuteron. 
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4.2 Trigger 

The trigger used for the recording of inelastic scattering events required signals from the 
hodoscopes HO, HI, and H2 and a sufficiently large energy deposition in the calorimeter. 
The efficiencies of the trigger detectors are extracted individually from special calibration 
triggers and combined to obtain the total trigger efficiency. It is assumed that inefficiencies 
in the electronic trigger logics are negligible. Trigger efficiencies are sensitive to, among 
others, misalignment effects, detector-voltage setting, and radiation damage, the latter 
especially in the HO hodoscope. This last effect is responsible for the reduced efficiencies 
seen at small scattering angles and for the differences between the top and bottom detector. 
Such differences are shown in Fig. |2| for data taken in the year 2000, for the kinematic 
binning used in the analysis. 
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Figure 2: Trigger efficiencies Etrigger for data taken in the year 2000 shown separately for the 
top and bottom spectrometer halves. The error bars represent only statistical uncertainties. The 
symbols refer to the bins shown in Fig. ^ 



4.3 Particle Identification 

The scattered lepton (positron or electron) is identified by a combination of the responses 
of the transition-radiation detector TRD, the pre-shower hodoscope H2, and the lead- 
glass calorimeter. Each of these elements used alone gives a high rejection of hadrons. A 
Cerenkov detector provides additional hadron identification. (A threshold Cerenkov was 
used for pion identification in 1996-97, and a ring-imaging Cerenkov detector was used 
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thereafter to identify pions, kaons, and protons.) The detector response of an individual 
PID element is determined by placing very restrictive constraints on the response of the 
remaining elements, thereby generating a clean sample of a given kind of particles with 
which the unit under study is calibrated. 

In combination, the array of detectors provides an average lepton identification with 
an efficiency Siep of about 98% and an average hadron contamination Chad below 1% over 
the full kinematic range of the HERMES acceptance. At low values of x, lepton efficiencies 
as low as 94% and hadron contaminations as high as 3% are reached. 

The efficiency for lepton identification and the fractional hadron contamination as a 
function of x for the various bins are presented in Fig. ^ for representative data taken 
in the year 2000. The figure shows that for smaller values of x (x < 0.1) a lower lepton 
identification efficiency appears correlated to a larger hadron contamination. 




Figure 3: Lepton identification efficiency £iep and hadron contamination Chad in the year 2000. 
The symbols refer to the binning shown in Fig. |l|. 



4.4 Charge-symmetric background 

The observed event sample is contaminated by background coming mostly from charge- 
symmetric processes, such as meson Dalitz decays (e. g. vr*^ — )• e+e~7) or photon conversions 
into e^e~ pairs. Since these positrons and electrons originate from secondary processes, 
they typically have lower momenta and are thus concentrated at high y. A correction for 
charge-symmetric background events Ncs is applied in each kinematic bin by counting with 
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negative weight leptons with a charge opposite to that of the beam particle. It is assumed 
that acceptance and inefficiencies are the same for background electrons and positrons, even 
though their spatial distributions after the magnet are quite different. The x dependence 
of CS, the ratio of charge-symmetric events to the total number of events in each kinematic 
bin, is shown for the six bins in Fig. The charge-symmetric background is negligible 
at large particle momenta, but reaches up to 12% at low particle momenta of about 6 GeV. 

4.5 Luminosity 

The integrated luminosity L per nucleon is calculated as flollows: 

L= I Cdt= {RlR - 2At ■ Rl ■ Rr) ■ Ctoe • CLumi ■ A6 ■ 



A 
Z 



(4.3) 



Here, Rl and Rr are the count rates in the left and right luminosity detector, respectively, 
Rlr is the coincidence rate measured within a time window of At = 40 ns, qj„e is the trigger 
livetime factor, Ciumi is the year-dependent luminosity factor, A6 is the time interval in 
which the luminosity rates were obtained, and A/Z is the ratio of the numbers of nucleons 
(A) and electrons (Z) in the target gas atoms. 
The term 2At -Rl-Rr in Eq. 
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corrects for accidental coincidences ac- 
cording to the statistical expectation and 
is of the order of 0.1-0.5%. The physics 
trigger livetime contribution cu^e is de- 
fined as the fraction of the physics events 
that are accepted by the data acquisi- 
tion system out of all events generating 
a physics trigger. This quantity is typi- 
cally above 90%. 

The data aquisition system of the 
luminosity detector worked independently 
of the physics triggers. It is assumed 
that the inefficiency of the luminosity 
event trigger was negligible. The lumi- 
nosity factor CLumi accounts for the ge- 
ometric acceptance of the luminosity de- 
tector, the beam position and the abso- 
lute M0ller and Bhabha cross sections. 
Its year dependence derives from the age- 
ing of the luminosity detector and dif- 
ferent running conditions, i.e., changes in beam charge and beam optics. The dependence 
of the coincidence rate from beam position and slope was measured in order to disentangle 
the dependence of the measured coincidence rate from beam orbits and geometrical ac- 
ceptance of the luminosity detector. The uncertainty in the measurement of the absolute 
luminosity is dominated by the uncertainty on the acceptance of the detector, which de- 
pends sensitively on the impact coordinates of the particle. The uncertainty on the latter 



Figure 4: Percentage of charge-symmetric back- 
ground, calculated from the ratio of the charge- 
symmetric events to the total events in each bin, for 
the 2000 deuterium data. The symbols refer to the 
binning shown in Fig. |l|. 
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is about 2.5 mm, which propagates into an uncertainty of about 7% on the integrated cross 
section and therefore on the luminosity. 



4.6 Instrumental smearing and radiative effects 

Instrumental smearing is due to intrinsic detector resolution and multiple scattering in the 
various detector elements of the particles emerging from the DIS process and identified as 
the scattered lepton. Radiative effects include vertex corrections to the QED hard scat- 
tering amplitude and radiation of one or more real photons by the incoming or outgoing 
lepton. Radiative effects and instrumental smearing both modify the Born kinematic con- 
ditions resulting in altered reconstructed kinematic variables. Migration probabilities for 
the relevant kinematic variables are determined from a Monte Carlo simulation and used 
to correct the measured distributions. 

The Born cross section for inelastic scattering on the proton is simulated according to 
the ALLM97 parameterization of F2 (2^ and the parameterization -R1990 for R, while 
that for the deuteron is derived from the same parameterizations in conjunction with the 
fit [||] to F^/F^ data from NMC, SLAC and BCDMS. Radiative effects are simulated with 
RADGEN ||3^. The electric and magnetic form factors of the proton and neutron, from 
which the elastic cross sections are derived, are taken from the fits in Refs. and p7| . 
When using the more recent parameterizations for the proton from Ref. |38|, the results 
are essentially the same. 

The probabilistic information about event migration can be summarized in a smearing 



matrix [39|, 



Si,J) = p2^^^^ (4.4) 

Here, n.{i,j) is the migration matrix representing the number of events originating from 
kinematic bin j at Born level and measured in bin i. It is extracted from a Monte Carlo sim- 
ulation with full track reconstruction that simulates the inelastic scattering cross section, 
QED radiative effects and instrumental smearing. Material outside the detector acceptance 
is excluded from this simulation for computational economy. The vector nBomij) contain- 
ing the number of events at Born level is obtained from a second Monte Carlo calculation 
that simulates only the (unradiated) inelastic cross section. An additional column j = 
is defined for events that migrate into the acceptance from outside. The smearing matrix 
S{i,j) has the property of being independent from the generated cross section within the 
acceptance. 

The inverted squared submatrix S'{i,j) = S{i,j > 0) relates the measured distribu- 
tions to the distributions at Born level: 

CTBornij) = '^S''^{j,i) X [aExpii) " S{i,0)aBorniO)]. (4.5) 

i 

The reconstruction of simulated tracks uses the same algorithm as for real data. 
Tracking-related inefficiencies are taken into account in the unfolding procedure, assuming 
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that coincident particles outside the acceptance do not significantly affect the efficiency 
and the simulation adequately models the physical processes in the tracking detectors. 

The Monte Carlo generated data sample was a factor 10 larger than the experimental 
data sample. The statistical uncertainties of the Monte Carlo data enter mainly via the 
simulated experimental count rates in the migration matrices. A multi-sampling numer- 
ical approach is used to propagate these statistical uncertainties through the unfolding 
algorithm. The statistical uncertainties of the inelastic scattering Born cross section com- 
ing from the experimental cross section and those originating from the finite statistical 
precision of the Monte Carlo are summed in quadrature to produce the total statistical 
uncertainty. 



4.7 Detection efficiency of specific radiative events 

Radiative corrections include cases where the incoming electron radiates a high-energy 
photon and then scatters elastically from the nucleon with negligible momentum transfer. 
The efficiency to detect such events is reduced due to the following effect. The radiated 
photon is emitted at small angles and has a large probability to hit the beam pipe, gen- 
erating an electromagnetic shower that saturates the wire chambers. This makes the data 
acquisition system skip the event as no tracking is possible. In order to compensate for this 
omission, the detection efficiency for elastic and quasi-elastic radiative events is estimated 
using a dedicated Monte Carlo simulation that includes a complete treatment of showers 
in material outside the geometrical acceptance. 

The resulting efficiencies Ee.m. are significantly less than 100% in the range 0.01 < 
X < 0.1. They show a dip at x ~ 0.02 where, in the case of the proton, they reach values 
as low as 80% while in the case of the deuteron they are as low as 90% (60%) for elastic 
(quasi-elastic) events. They are applied to the background term S{i,0)aBorn{0) in order 
to not over-correct for radiative processes that are not observed in the spectrometer. More 



details can be found in Refs. |39|, 40 1. 



4.8 Misalignment effects 

Imperfect alignment of the two spectrometer halves and the beam with respect to their 
ideal positions is studied in order to estimate the impact on the measured cross sections and 
structure functions. Misalignment effects cannot be corrected for in the unfolding because 
they are not of a stochastic nature. Rather, they are studied in a Monte Carlo simulation, 
and the fractional change of the Born cross section in each kinematic bin is obtained from 
the ratio of unfolded cross sections when using a MC with an aligned geometry and another 
with a misaligned geometry. These fractional changes are used to rescale the experimental 
Born cross sections on a bin-by-bin basis. Misalignment effects are most significant for 
small scattering angles and high particle momenta, i.e., at small in each bin of x. The 
correction reaches values as high as 19% in the lowest bin and decreases to about 3% 
in the highest bin. 
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5. Systematic Uncertainties 



5.1 Inclusive inelastic scattering cross sections 

Particle identification. 

Correlations between PID detectors as described in Sect. 4.2 cannot be completely avoided 



They may bias the correction for particle identification. These effects are covered by the 



assignment of a conservative PID uncertainty of the full size of the correction (see Eq. (4.1)). 
Hadrons are predominantly produced at small momenta. Thus particle misidentification 
occurs more likely at high y, i.e., towards higher in each bin of x. Nevertheless, 
the uncertainty due to particle identification, 6pjd, is always smaller than 3%, because 
contaminations somewhat compensate inefficiencies. 



Instrumental smearing and radiative effects. 

In the unfolding procedure an uncertainty can arise from uncertainties in the formalism to 
calculate radiative effects and in the model used for the cross section outside the accep- 
tance. The latter affect our results through the radiative tail. The uncertainty, Smodeh was 
estimated by varying the input elastic and inelastic cross sections within their uncertain- 
ties and found to be below 2%, except for a few bins, where it went up to 4.3% (3.1%) 
at maximum for the proton (deuteron) case. This is negligible compared to the overall 
normalization uncertainty of our data of about 7% (see below). 

Misalignment. 

In each bin, half of the deviations of yields obtained in a Monte Carlo simulation with 
estimated geometric misalignments from the yields obtained in a Monte Carlo simulation 
with aligned (ideal) geometry serve as an estimate of the systematic uncertainty due to 
misalignment. The uncertainty due to misalignment, 6mis., reaches values of up to 5.4%. 
However, the majority of the data points has an uncertainty due to misalignment that is 
smaller than 3%. 



Dependence on misalignment of the efficiency Se.m. for elastic and quasi-elastic radiative 
events. 

The efficiency £e.m. and its dependence on misalignment were studied in Monte Carlo 
simulations. The assignment of a corresponding systematic uncertainty is accomplished 
by applying the values of £e.m. extracted from a Monte Carlo simulation with aligned and 
misaligned geometry to the high-multiplicity radiative events included in the background 
term S{i, 0)crBorn{0). The difference of the unfolded inelastic scattering Born cross sections 
obtained for these efficiencies is assigned as a systematic uncertainty, 6rad.i due to these 
radiative corrections. 



Overall normalization uncertainty. 

The normalization uncertainties of the absolute cross sections and the structure functions 
are dominated by the uncertainty of the year-dependent luminosity constant Ciumi in 
Eq. ( [4.3[ ) . The uncertainties of the luminosity constants weighted with the sizes of the data 



- 13 - 



sets result in an overall normalization uncertainty of 7.6% for the data taken on a hydrogen 
target and 7.5% for the data taken on a deuterium target. 

5.2 Inclusive inelastic scattering cross-section ratio a'^ /g'p 

The cross-section ratio can be determined with higher precision than the cross sections 
themselves due to the cancellation of the misalignment uncertainty, the PID uncertainty 
and, to a large extent, the overall normalization uncertainty. The remaining overall nor- 
malization uncertainy of 1.4% is attributed to variations of the beam conditions between 
data sets. 

The efficiencies £e.m. for proton and deuteron are different and therefore do not 
cancel in the proton-to-deuteron cross-section ratio. The uncertainty dj-ad of the cross- 
section ratio is obtained by propagating, for proton and deuteron cross sections, the un- 
certainties of efficiencies for high- multiplicity radiative events due to misalignment. It is 
found to be less than 2.5% in every kinematic bin. 



6. Discussion of the Results 

The kinematic conditions of the HERMES inclusive lepton-nucleon scattering cross sections 
presented here overlap those of existing data over a large kinematic range. New information 
is provided in the region with '^1 GeV^ and 15GeV^ < W"^ < 45GeV^, corresponding 
to 0.006 <x < 0.04. 

6.1 Structure functions and 

The differential cross sections dPa^''^ / dxdQ'^ for inelastic scattering on the proton and 
deuteron as well as the corresponding structure functions and Ff are listed in Tabs. ^ 
and |6| together with the statistical and systematic uncertainties. The statistical uncertain- 
ties of the HERMES data range between 0.4% and 3.0%. Almost 80% of all data points 
have a statistical uncertainty smaller than 1%. The overall normalization uncertainty of 
7.6% (7.5%) for the inelastic scattering cross section on the proton (deuteron) and the 
contribution from misalignment are the dominating systematic uncertainties. 

The differential cross sections are shown in Figs. || and|6| as a function of in bins of x. 
The structure functions are shown in Figs. ^ and |8|, together with the available world data 
from fixed target (E665 |l8|, BCDMS NMC 0, SLAG JLAB |0|, ^, |2|, ||) 



and collider experiments (HI and ZEUS). The data are overlaid with new fits to world data, 
including the data presented here, of inclusive proton (GDll-P) and deuteron (GDll-D) 
cross sections. These functions are described in section |6.3| . 

In the region x > 0.07 and > \ GeV^, HERMES data are in good agreement with 
existing data from SLAG and NMG. The HERMES measurement provides also data in a 
previously uncovered kinematic region between JLAB data on the one hand and NMG, 
BGDMS, E665 and the collider experiments on the other. This can be clearly seen in 
Figs. and |[ 
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For virtual-photon data (Q^ > 0), the photon- nucleon cross section (t^'^-j, = a^jf + a^' 
can be derived from the structure function Fl'''^ using Eq. ( ^.41) : 



<^ = 4vr^aen. ^^,,; Ff'^ . (6.1) 



The W'^ dependence of the resulting proton and deuteron cross sections o'^'^y is shown in 
Figs. ^ and |l^ for HERMES results on inelastic scattering in comparison with world data 
and real-photon cross sections. HERMES data fill a gap for in the range 0.3-1. GeV^ 
and W"^ in the range 20-50 GeV^. 

6.2 Cross-section ratio a'^/a^ 

From the measured cross-section ratio, a'^ /a^ , the ratio of the deuteron and proton struc- 
ture functions, F2 /F2, can be extracted. In the DIS regime F2/F2 is related to the ratio of 
the down and up quark distribution functions and imposes strong constraints on the x de- 
pendence of the flavor composition of the nucleon. From the combination of F2 and F2 /F2 
the non-singlet structure function F2 — F2 and the Gottfried sum Sq = / {F2 — Fll^dxjx 
can be determined. For the latter, it was shown that the simple quark-model expectation 



of 1/3 is not reached |42, 43 1 and hence the light quark sea is not flavor-symmetric. We do 
not elaborate on these aspects here, as our data add only limited additional information 
in the DIS region. 

The cross-section ratio cr'^/o"*' is determined year-by-year and then averaged. Thereby 
the normalization uncertainty is reduced and the effects of PID efficiencies and contami- 
nations cancel. Because the deuteron and the proton data were partly taken with different 
beam conditions and resulting possible differences in absolute normalization, the ratio of 
the -F2-values given in Tables |5| and ^ may differ slightly from the ratio values presented 
here. The results for ct'^/cj^ are listed in Tab. ^. The values are shown as a function of 



in bins of x in Fig. 11 together with world data. 

6.3 Fits to world data to the cross section o^lt^rp. 

Fits of the inclusive inelastic scattering cross sections on the proton and deuteron are 
performed. The ALLM functional form selected for the fits is described in Refs. Q and p5|| . 
It is a 23-parameter Regge-motivated model of o^^^j. = a^jf + a^'^ ^ valid for W'^ > 4GeV^, 
i.e., above the resonance region, and any including the real photon point (Q^ = 0). 

These new fits to the photon-nucleon cross sections (tl+t on proton (GDll-P) and 
deuteron (GDll-D) reflect the recent world knowlege on the cross sections. In order to 
present a self-consistent set, the values of F2 to be used in the fit are calculated from the 
measured cross sections using for all data the parametrization -Riggs. In cases where the 
measured cross sections were not given, the latter are reconstructed from the published 
values for F2 using the values of R used in calculating these F2. Both fits GDll-P and 
GDll-D are performed in the same manner. 

The fit to proton data includes 2441 data points from the SLAG experiments E49a, 



E49b, E61, E87, E89a, E89b p|, NMG Q, the combined HERA data [Q, E665 []l8|. 



BGDMS [|15|, JLAB |2g, ||, Hi, the HERMES data presented here, and real photon 
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data |45|. In the case of HERA, only data points with < 1000 GeV^ were used in order 
to minimize contributions from weak processes to the cross section. The fit to deuteron 
data includes 2497 data points from the SLAG experiments E49a, E49b, E61, E87, E89a, 
E89b, E139, E140 g9|, NMC 0; E665 BCDMS @; JLAB [|l],|2|,||; the HERMES 
data presented in this paper, and real photon data |45]. Tab. |3| lists in more detail the data 
sets, together with their (x) and (Q^) ranges. 



Experiment 



# of 


{x) range 


{Q ) range 


target 


points 




[GeV J 




456 


0.621x10 - 0.65 


0.045 - 1000 


P 


95 


0.13x10 - 0.65 


90 - 1000 


P 


91 


0.89x10 - 0.39 


0.23 - 64 


p, d 


73 


0.78x10 ^ - 0.465 


0.8 - 8.75 


p, d 


65 


0.87x10 ^ - 0.475 


1.18 - 14.4 


p, d 


75 


0.35x10 ^ - 0.477 


0.83 - 34.8 


p, d 


79 


0.37x10 ^ - 0.479 


1.27 - 62.3 


p, d 


58 


0.07 - 0.75 


7.5 - 75 


P 


62 


0.07 - 0.75 


8.75 - 99 


p, d 


57 


0.07 - 0.75 


8.75 - 137 


P 


56 


0.07 - 0.75 


17 - 137.5 


d 


52 


0.1 - 0.75 


32.5 - 230 


p, d 


98 


0.067 - 0.71 


0.59 - 7.9 


p, d 


187 


0.08 - 0.854 


0.663 - 20 


P 


174 


0.08 - 0.854 


0.663 - 20 


d 


25 


0.065 - 0.326 


0.58 - 1.68 


P 


24 


0.065 - 0.326 


0.58 - 1.68 


d 


94 


0.313 - 0.832 


3.96 - 19.7 


p, d 


72 


0.31 - 0.9 


3.6 - 30.3 


P 


66 


0.33 - 0.9 


3.86 - 30.2 


d 


98 


0.063 - 0.809 


0.887 - 18.4 


P 


81 


0.063 - 0.807 


0.887 - 18.4 


d 


18 


0.089 - 0.7 


2 - 15 


d 


38 


0.2 - 0.5 


1 - 10 


d 


50 


0.483 - 0.641 


3.585 - 5.67 


p, d 


272 


0.23 - 0.509 


1.325 - 3.275 


P 


1018 


0.122 - 0.617 


0.475 - 5.125 


d 


5 


0.25 - 0.8 


0.15 - 1.045 


p, d 


50 


0.9x10^2 - 0.25 


0.034 - 1.761 


p, d 


81 


0.69x10-2 - 0.664 


0.291 -12.78 


p, d 


196 






P 


174 






d 



HERA (positron beam) 
HERA (electron beam) 
E665 @ 

NMC - 90 GeV @ 
NMC - 120 GeV |7| 



NMC - 200 GeV |T7| 
NMC - 280 GeV |7| 
BCDMS - 100 GeV |5| 



BCDMS - 120 GeV |5| 
BCDMS - 200 GeV |5| 
BCDMS - 200 GeV |5| 
BCDMS - 280 GeV p| 
SLAC E49a |l|] 
SLAG E49b p| 
SLAC E49b 



SLAC E61 191 



SLAC E61 
SLAC E87 m 
SLAC E89a 
SLAC E89a ^ 
SLAC E89b ^\ 
SLAC E89b p| 
SLAC E139 H 
SLAC E140 H 
JLAB EOO-116 



JLAB CLAS H 



JLAB CLAS |2| 

JLAB HALL C (Rosenbluth) [|| 

JLAB HALL C (Model Dep.) H 

HERMES (this analysis) 

real photon ^5| 



real photon |45| 



Table 3: Data sets used in the GDll fits. 
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The fits are based on the minimization of the value of defined as 

i,k y^i^k + '^i,k i u + °kVk) k 

^ sr^ [D^AW\ Q') - T{p, W\ Q-) • (1 - 6^1^^)? ^ ^,2 .q^. 

i,k ^i,k +^i,k k 

where Dj^^ it u**^* '^^l^k^ ^^^^ the values of gl+t for data point i within the data set k, 5k is 
the normalization uncertainty in data set k quoted by the experiment, Vk is a normalization 
parameter, T(p,VF^,Q-^) is the 23-parameter ALLM functional form, p is the vector of 
functional parameters and ly is the vector of normalization parameters. The definition of 
takes into account point-by-point statistical and systematic uncertainties. We agree 
that ideally one should treat the statistical and systematic uncertainties separately, as the 
latter in most cases are to some degree correlated. However, almost all published data 
sets just give statistical, total, and normalization uncertainties. For that reason we have 
used in the fit the total uncertainties per point, plus the overall normalization. Clearly the 
value of is influenced by this, but in practice the influence on the resultant parameters 
seems minor, see for example Ref. |46|. Additionally, for each data set k a normalization 



parameter u^. is used to describe the normalization of all data points belonging to the data 
set. Since the normalization parameters are assumed to normalize data within the known 
normalization uncertainties, they are scaled with these uncertainties and taken into account 
by a penalty term v'^ to control their variation according to their standard deviations. 

At each iteration of the x^-minimization, the normalization parameters are analyt- 
ically determined: 

Ei hTi,k{Ti,k - Di^k)/o-lk 
2li Ti^k^k/'^^k + 1 

where cr^^ = crf^^'^ + (jf^k^ . This equation is obtained by requiring dy^ jdvk = in 
the context of the approximation for in Eq. (|6.2|). The analytical determination of the 
parameters Vk reduces the number of free parameters in the x'^-niinimization. The separate 
extraction of all normalization parameters is possible only because the parameters Vk are 
uncorrelated. At each step Eq. (|6.3| ) is substituted into the x^-function, Eq. (|6.2| ), and 
thus the final fit result is obtained by minimizing x^ only with respect to the functional 
parameters. Similar methods are used by others in order to save computing time 48]. 



The fits are performed using the CERNLIB package MINUIT Q and the MINUIT exten- 
sion |5C]. The resulting fit parameters are listed in Tab. ^ together with the parameters 
obtained in the previous fits ALLM97 j25] and GD07-P performed on proton data. 
The difference between GD07-P and GDll-P is the presence of HERMES, JLAB and the 
combined HERA data in the most recent fit. In Ref. [ETI the uncertainties were calculated 



with an UP value |4£] equal to 24.7. That identical fit was redone here with the same value 
of UP=1 as the new fits in order to provide an appropriate comparison of its uncertain- 
ties. The ALLM97 fit provided no uncertainties. The uncertainties given in the table for 
GD07-P, GDll-P, and GDll-D correspond to the diagonal elements of the full covariance 
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matrix, which must be used to calculate uncertainties in F2 or cross sections. The / per 
degree of freedom for the GDll-P is 0.89 (it is 0.92 for GD07-P), while for the deuteron fit 
GDll-D it is 0.68. The fitted relative normalization of HERMES proton (deuteron) data 
is +1.5% (-2.2%), well within the quoted uncertainty of 7.5% (7.6%). Tab. ^ shows the 
fitted relative normalizations for all the data sets used in the fit. 

6.4 Fits to world data for the cross-section ratio a'^ /a^. 

The world data on the cross-section ratio a'^ /a^ is fitted using the functional form 

va'^/aP{x, Q2) = A{x) + B{x) In (6.4) 

with 

A{x) = Pl+ P2X + P3X^ + pax^ + P5X'^, 

B{x) = Pa + p-jx + psx'^ , (6.5) 

previously used by NMC in Ref . [ 34 1 . It includes eight free parameters pi, ...,ps. The full 



data set consists of 260 data points from NMC ||5^, 621 data points from SLAG |1£], 49 data 



points from JLAB [23|, 81 data points from this analysis, 159 data points from BGDMS, 



and 53 data points from EMG. The latter two data sets are taken from Ref. |19]. The 
relative normalizations are also fitted using the same method as described in the previous 
section. The final results are listed in Tab. ^. 

The relative normalizations were all found to be well within the normalization uncer- 
tainties quoted by the experiments. In the case of HERMES data, it is found to be 0.996. 
Fig. |ll| shows the world data of a'^/a^ overlapped with the fits described in this section 
and with the ratio of to i^f derived from the GDll fits. While the ratio deviates from 
unity at large x, both GDll and NMC-like fits describe the data well but start to deviate 
from each other outside the range where the ratio has been determined directly. 



7. Summary 

This high-statistics measurement of the inelastic scattering cross section on the proton and 
the deuteron by HERMES provides data in the ranges 0.006 < x < 0.9 and 0.1 GeV^ < 

< 20 GeV^ and contributes data in a kinematic region previously not covered, ^ 1 
GeV^ and ISGeV^ < W"^ < AbGeY^, corresponding to 0.006 <x < 0.04. In the region of 
overlap, HERMES data agree with the world data. The new results are fitted in conjunction 
with the other data, including those at the photon point, to give a new parameterization 
of the structure functions (^2^) for = to 30000 GeV^ (230 GeV^) and x-values 
0.6 X 10"^ (0.89x10-^) to 0.9. The fitted relative normalizations of HERMES data are 
well within the quoted normalization uncertainties. 

The proton and deuteron cross-section ratio was also determined. This quantity has 
a much smaller point-to-point uncertainty because several systematic uncertainties cancel 
in the ratio and its normalization uncertainty is also much smaller. The HERMES data 
show excellent agreement with the other data and the fitted normalization is well within 
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the quoted normalization uncertainty. The cross-section ratio from HERMES is also fitted 
in conjunction with existing data from other experiments. The dependence can be well 
described by a simple (x-dependent) In term. 
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parameter 


value 


uncertainty 


V\ 


0.986 


0.0021 


P2 


-0.684 


0.034 


P3 


1.54 


0.18 


PA 


-2.93 


0.36 


Pb 


1.96 


0.22 


P6 


0.0053 


0.0011 


P7 


-0.0974 


0.0087 


IM 


0.125 


0.014 


normalization 


value 




HERMES 


0.996 




NMC 


0.999 




BCDMS 


1.010 




,- SLAG 


1.003 




e JLAB 


1.000 




EMC 


0.995 





Table 4: Results of fit to the cross-section ratio 

a'^/o'P. The meaning of the parameters pi,...ps 
and relative normalizations is explained in the 
text. 
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Figure 5: Inelastic proton differential DIS cross section cPaP/dxdQ'^, multiplied by a factor 
for the purpose of illustration, in the kinematic range 0.008 < (x) < 0.639 and 0.2 GeV^ < 
(Q^) < 20 GcV^. The values of (Pa^ /dxcLQ"^ ■ Q'^ are scaled by powers of 1.6. The results are 
overlaid with the phenomenological parameterization GDll-P (central curves) and its uncertainty 
(outer curves). Details on the fits can be found in Sect. The error bars represent the total 
uncertainties calculated as the sum in quadrature of all statistical and systematic uncertainties 
including normalization, and are smaller than the symbols. 
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Figure 6: Inelastic deuteron differential DIS cross section <P a'^ / dxdQ"^ , multiplied by a factor 
Q"^ for the purpose of illustration, in the kinematic range 0.008 < {x) < 0.639 and 0.2 GeV^ < 
(Q^) < 20 GcV^. The values of d'^a''- /dxdQ'^ ■ Q"' are scaled by powers of 1.6. The results are 
overlaid with the phenomenological parameterization GDll-D (central curves) and its uncertainty 
(outer curves). Details on the fits can be found in Sect. The error bars represent the total 
uncertainties calculated as the sum in quadrature of all statistical and systematic uncertainties 
including normalization, and are smaller than the symbols. 
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Figure 7: HERMES data for F2 together with world data in the kinematic range 0.008 < (x) < 
0.679 and 0.02 GeV^ < (Q^) < 20 GeV^. The results are overlaid with the phenomenological 
parameterization GDll-P (black solid central curve) and its uncertainty (outer curves) obtained as 
described in Sect. 6.3. A bin-centering correction is applied to the data in order to match the central 



values of the x bins. The values of F2 are scaled by powers of 1.6. Inner error bars are statistical 
uncertainties, while outer error bars are total uncertainties calculated as the sum in quadrature of 
all statistical and systematic uncertainties including normalization. 
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Figure 8: HERMES data for F2 together with world data in the kinematic range 0.008 < (x) < 
0.679 and 0.02 GeV^ < (Q^) < 20 GeV^ The results are overlaid with the phenomenological 
parameterization GDll-D (black solid central curve) and its uncertainty (outer curves) obtained as 
described in Sect. 6.3. A bin-centering correction is applied to the data in order to match the central 



values of the x bins. The values of F2 are scaled by powers of 1.6. Inner error bars are statistical 
uncertainties, while outer error bars are total uncertainties calculated as the sum in quadrature of 
all statistical and systematic uncertainties including normalization. 
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Figure 9: HERMES data for the photon-proton cross section function of , together 

with world data and the results from the GDll-P fit (central curves) and its uncertainties (outer 
curves), in bins of Q^. The data points denoted 'real photon' are for photoproduction. Inner error 
bars are statistical uncertainties, while outer error bars are total uncertainties calculated as the 
sum in quadrature of all statistical and systematic uncertainties including normalization. 
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Figure 10: HERMES data for the photon-dcutcron cross section function of W , 

together with world data and the resuUs from the GDll-D fit (central curves) and its uncertainties 
(outer curves), in bins of Q^. The data points denoted 'real photon' are for photoproduction. Inner 
error bars are statistical uncertainties, while outer error bars arc total uncertainties calculated as 
the sum in quadrature of all statistical and systematic uncertainties including normalization. 
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Figure 11: The dependence of a'^/a" for data from SLAC, HERMES, EMC, NMC, JLAB 
and BCDMS in bins of (a;). The error bars on the data points represent tot al uncertainties. The 
data points are overlaid with the results from the fit described in section 3.2 (blue solid line), with 
the unity line (red dashed) and with the ratio of the GDll-D over GDll-P fits (green dash-dotted 
line) . 
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A. Tables 
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e 
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r 
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1 


A 


0.0069 
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0.323x10^ 
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1.60 


0.69 


4.30 


5.36 
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0.0102 


0.404 
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1.35 
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3.95 
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0.514 
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0.0186 
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Table 5: Results on the differential Born cross section f^^Qi and . The statistical uncertainty 
5stat. and the systematic uncertainties Spiu (particle identification), Smodei (model dependence 
outside the acceptance), Smis. (misalignment), and 6rad. (Bethe-Heitler efficiencies) are given in 
percent. Corresponding x bin numbers and bin numbers and the average values (x) and (Q^) 
are listed in the first four columns. The overall normalization uncertainty is 7.6%. The structure 
function F2 is derived using the parameterization R = Rig^s ]3^ - 
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0.04 


0.27 


1.91 


0.08 


15 


D 


0.2106 


3.709 


0.219x10^ 


0.343 


0.68 


0.05 


0.22 


1.80 


0.09 


16 


D 


0.2726 


4.527 


0.106x10^ 


0.315 


0.84 


0.05 


0.19 


2.14 


0.02 


1 7 


pi 




P, 7Q7 


U.ZDo X lU 


9/19 


u.o / 




n 1 

U. iO 


1 .Do 


U.Ui 


18 


D 


0.5038 


11.344 


0.350 


0.132 


0.98 


0.07 


0.32 


1.61 


0.01 


9 


E 


0.0597 


1.906 


0.210x10^ 


0.333 


0.93 


0.30 


0.87 


0.78 


0.17 


10 


E 


0.0726 


2.276 


0.128x10^ 


0.346 


0.93 


0.22 


0.62 


1.53 


0.16 


11 


E 


0.0885 


2.587 


0.863x102 


0.351 


0.70 


0.09 


0.47 


1.33 


0.03 


12 


E 


0.1083 


3.071 


0.519x10^ 


0.355 


0.58 


0.03 


0.37 


1.49 


0.01 



Table 5: - continued. 
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X bin 


(c^ Din 


{X) 


/n2\ 
\v / 


dx dQ"^ 




Ostat. 


OPID 


A , , 

"model 


A 

'Jmis. 


A 

'-'rad. 








[GeV2] 


[nb/GeV2] 




[%] 


[%] 


[%] 


[%] 


[%] 


13 


E 


0.1335 


3.538 


0.331x102 


0.356 


0.82 


0.02 


0.29 


0.91 


0.01 


14 


E 


0.1664 


4.051 


0.210x10^ 


0.354 


0.56 


0.04 


0.23 


1.32 


0.04 


15 


E 


0.2107 


4.702 


0.122x102 


0.336 


0.58 


0.05 


0.18 


1.02 


0.05 


16 


E 


0.2726 


5.685 


0.603x10 


0.304 


0.66 


0.04 


0.14 


2.03 


0.03 


17 


E 


0.3679 


9.154 


0.124x10 


0.234 


0.68 


0.05 


0.27 


1.40 


0.01 


9 


F 


0.0603 


2.319 


0.124x103 


0.330 


0.76 


2.40 


1.35 


0.77 


0.13 


10 


F 


0.0734 


2.806 


0.726x102 


0.339 


0.69 


2.17 


1.01 


0.48 


0.15 


11 


F 


0.0894 


3.341 


0.444x102 


0.349 


0.68 


1.58 


0.77 


0.66 


0.11 


12 


F 


0.1094 


3.946 


0.267x102 


0.346 


0.62 


0.96 


0.60 


1.03 


0.08 


13 


F 


0.1347 


4.618 


0.164x102 


0.346 


0.58 


0.41 


0.46 


0.56 


0.03 


14 


F 


0.1679 


5.344 


0.102x102 


0.342 


0.58 


0.11 


0.34 


1.45 


0.04 


15 


F 


0.2128 


6.266 


0.592x10 


0.330 


0.63 


0.02 


0.27 


2.05 


0.02 


16 


F 


0.2759 


7.501 


0.299x10 


0.296 


0.57 


0.04 


0.24 


0.85 


0.01 


Table 5: - continued. 


j; bin 


Q2 bin 


(x) 


m 


^2(7^ 

dx dQ'^ 


Fi 


Sstat. 


SpiD 


Srnodel 


Smis. 


Srad. 








[GeV2] 


[nb/GeV2] 




[%] 


[%] 


[%] 


[%] 


[%] 


1 


A 


0.0069 


0.291 


0.327x10^ 


0.164 


1.41 


0.69 


3.14 


5.36 


0.36 


2 


A 


0.0102 


0.404 


0.138x10^ 


0.193 


1.30 


0.87 


2.15 


3.95 


0.47 


3 


A 


0.0142 


0.514 


0.724x10^ 


0.217 


0.84 


0.61 


1.70 


4.36 


0.40 


4 


A 


0.0186 


0.626 


0.435x10^ 


0.242 


0.95 


0.36 


1.53 


3.89 


0.30 


5 


A 


0.0244 


0.725 


0.304x10^ 


0.274 


0.81 


0.10 


1.43 


4.66 


0.27 


6 


A 


0.0325 


0.872 


0.178x10"^ 


0.291 


0.87 


0.01 


1.38 


3.78 


0.22 


7 


A 


0.0398 


0.903 


0.155x10^ 


0.306 


1.32 


0.05 


1.50 


3.92 


0.10 


8 


A 


0.0486 


0.929 


0.136x10^ 


0.323 


2.26 


0.07 


2.63 


6.43 


0.29 


9 


A 


0.0595 


1.102 


0.803 xlO^ 


0.323 


2.60 


0.05 


1.33 


3.17 


0.70 


10 


A 


0.0724 


1.127 


0.722x10^ 


0.348 


2.62 


0.06 


1.14 


5.49 


0.45 


11 


A 


0.0882 


1.167 


0.564x10^ 


0.338 


2.53 


0.07 


0.84 


3.45 


1.18 


12 


A 


0.1078 


1.205 


0.460 xlO^ 


0.345 


3.66 


0.07 


0.75 


5.93 


1.18 


13 


A 


0.1330 


1.263 


0.335x10^ 


0.329 


3.06 


0.07 


0.78 


2.56 


0.94 


14 


A 


0.1657 


1.337 


0.265x10^ 


0.354 


2.65 


0.06 


1.02 


7.52 


0.90 


15 


A 


0.2093 


1.447 


0.155x103 


0.298 


2.21 


0.06 


1.87 


2.99 


0.52 



Table 6: Results on the differential Born cross section f^^Q'i ^'^'^ ^2- The statistical uncertainty 
Sstat. and the systematic uncertainties Spijj (particle identification), S„iodei (model dependence 
outside the acceptance), Smis. (misalignment), and Srad. (Bethe-Heitler efficiencies), are given in 
percent. Corresponding x bin numbers and bin numbers and the average values (x) and (Q^) 
are listed in the first four columns. The overall normalization uncertainty is 7.5%. The structure 
function F2 is derived using the parameterization R = Riggs ||3^ • 
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X bin 


Q2 bin 


(x) 


[GeV2] 


dx dQ'^ 
[nb/GeV2] 




5stat. 
[%] 


^PID 

[%] 


^model 
[%] 


[%] 


^rad. 
[%] 


16 


A 


0.2642 


1.749 


0.845 xlO^ 


0.298 


1.25 


0.05 


2.65 


3.29 


0.56 


17 


A 


0.3477 


2.672 


0.196x10^ 


0.217 


1.37 


0.05 


1.99 


7.62 


0.21 


18 


A 


0.4605 


4.224 


0.377x10 


0.142 


1.13 


0.04 


1.70 


4.10 


0.11 


19 


A 


0.6346 


8.688 


0.198 


0.047 


2.80 


0.08 


1.25 


1.17 


0.02 


1 


B 


0.0085 


0.354 


0.203x10^ 


0.189 


1.29 


0.93 


2.43 


3.75 


0.81 


2 


B 


0.0117 


0.477 


0.979x10^ 


0.226 


0.77 


1.22 


2.05 


2.07 


0.69 


3 


B 


0.0152 


0.612 


0.488x10^ 


0.241 


0.61 


1.50 


1.90 


1.03 


0.36 


4 


B 


0.0193 


0.770 


0.260x10^ 


0.258 


0.64 


1.80 


1.79 


0.79 


0.40 


5 


B 


0.0250 


0.896 


0.169x10"^ 


0.272 


0.52 


1.07 


1.51 


0.56 


0.12 


6 


B 


0.0327 


1.144 


0.865 xlO^ 


0.291 


0.68 


0.92 


1.47 


0.64 


0.30 


7 


B 


0.0401 


1.074 


0.100x10^ 


0.308 


1.06 


0.01 


1.28 


2.61 


0.17 


8 


B 


0.0490 


1.123 


0.843x10^ 


0.318 


1.05 


0.03 


1.22 


2.47 


0.24 


9 


B 


0.0596 


1.281 


0.582 xlO^ 


0.337 


2.05 


0.06 


0.81 


2.26 


0.36 


10 


B 


0.0725 


1.359 


0.439x10^ 


0.329 


1.85 


0.06 


0.58 


3.96 


0.36 


11 


B 


0.0884 


1.399 


0.383x10^ 


0.349 


1.96 


0.06 


0.51 


3.58 


0.05 


12 


B 


0.1082 


1.554 


0.260x10^ 


0.348 


1.22 


0.06 


0.47 


4.34 


0.19 


13 


B 


0.1333 


1.632 


0.195x10^ 


0.339 


1.85 


0.06 


0.49 


1.84 


0.34 


14 


B 


0.1661 


1.811 


0.132x10^ 


0.341 


1.01 


0.05 


0.50 


4.10 


0.05 


15 


B 


0.2104 


2.050 


0.763x10^ 


0.314 


1.53 


0.05 


0.58 


2.80 


0.09 


16 


B 


0.2722 


2.471 


0.374x10^ 


0.285 


1.02 


0.05 


0.78 


3.85 


0.18 


17 


B 


0.3620 


3.790 


0.893x10 


0.218 


0.86 


0.04 


0.57 


2.21 


0.04 


18 


B 


0.4900 


6.026 


0.139x10 


0.120 


0.65 


0.04 


0.62 


0.94 


0.01 


19 


B 


0.6641 


12.778 


0.605x10-1 


0.036 


2.32 


0.05 


0.66 


0.86 


0.01 


5 


C 


0.0270 


1.092 


0.103x10"^ 


0.292 


0.68 


2.39 


1.81 


1.06 


0.25 


6 


C 


0.0338 


1.391 


0.535 xlO^ 


0.309 


0.99 


2.98 


1.93 


0.28 


0.43 


7 


c 


0.0401 


1.293 


0.621 xlO^ 


0.309 


0.80 


0.38 


1.35 


2.71 


0.15 


8 


c 


0.0490 


1.460 


0.430x10^ 


0.316 


0.70 


0.17 


1.05 


1.29 


0.12 


9 


c 


0.0596 


1.466 


0.400x10^ 


0.323 


1.09 


0.01 


0.69 


1.73 


0.24 


10 


c 


0.0726 


1.627 


0.291 xlO^ 


0.337 


1.10 


0.04 


0.48 


2.05 


0.07 


11 


c 


0.0885 


1.691 


0.236x10^ 


0.337 


1.11 


0.05 


0.41 


1.43 


0.02 


12 


c 


0.1082 


2.013 


0.140x10^ 


0.341 


0.96 


0.05 


0.37 


2.70 


0.08 


13 


c 


0.1334 


2.158 


0.102x10^ 


0.335 


1.09 


0.05 


0.33 


2.17 


0.05 


14 


c 


0.1662 


2.458 


0.636x10^ 


0.328 


1.22 


0.05 


0.30 


0.98 


0.07 


15 


c 


0.2105 


2.851 


0.361 xlO^ 


0.310 


0.75 


0.04 


0.29 


2.40 


0.07 


16 


c 


0.2724 


3.479 


0.167x10^ 


0.272 


0.67 


0.04 


0.28 


1.74 


0.03 


17 


c 


0.3622 


5.125 


0.428x10 


0.205 


0.62 


0.04 


0.21 


1.84 


0.04 


18 


c 


0.4957 


8.280 


0.630 


0.113 


0.70 


0.05 


0.22 


2.23 


0.02 



Table 6: - continued. 
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X bin 


Q2 bin 


(x) 


iQ') 

[GeV2] 


dx dQ"^ 
[nb/GeV2] 




5stat. 
[%] 


5piD 

[%] 


^model 
[%] 


[%] 


[%] 


7 


D 


0.0408 


1.615 


0.343 xlO^ 


0.311 


0.68 


2.58 


1.76 


0.53 


0.34 


8 


D 


0.0496 


1.927 


0.208 xlO^ 


0.320 


0.60 


2.49 


1.49 


0.68 


0.10 


9 


D 


0.0596 


1.675 


0.290x10^ 


0.328 


1.21 


0.04 


0.69 


1.88 


0.06 


10 


D 


0.0726 


1.928 


0.192x10^ 


0.339 


1.46 


0.01 


0.46 


2.13 


0.01 


11 


D 


0.0885 


2.080 


0.144x103 


0.339 


0.89 


0.03 


0.37 


1.21 


0.07 


12 


D 


0.1083 


2.490 


0.824x102 


0.336 


1.01 


0.04 


0.31 


1.50 


0.18 


13 


D 


0.1334 


2.795 


0.570x102 


0.345 


0.83 


0.04 


0.26 


2.67 


0.11 


14 


D 


0.1663 


3.195 


0.343x102 


0.326 


0.90 


0.04 


0.22 


1.91 


0.01 


15 


D 


0.2106 


3.709 


0.195x102 


0.305 


0.69 


0.05 


0.19 


1.80 


0.04 


16 


D 


0.2726 


4.527 


0.899x10 


0.266 


0.88 


0.05 


0.17 


2.14 


0.04 


17 


D 


0.3624 


6.797 


0.215x10 


0.197 


0.58 


0.04 


0.12 


1.63 


0.01 


18 


D 


0.5038 


11.344 


0.272 


0.103 


1.03 


0.07 


0.31 


1.61 


0.02 


9 


E 


0.0597 


1.906 


0.202x10^ 


0.321 


0.90 


0.30 


0.74 


0.78 


0.08 


10 


E 


0.0726 


2.276 


0.124x10^ 


0.334 


0.92 


0.22 


0.53 


1.53 


0.13 


11 


E 


0.0885 


2.587 


0.826x102 


0.335 


0.68 


0.09 


0.40 


1.33 


0.11 


12 


E 


0.1083 


3.071 


0.494x102 


0.338 


0.55 


0.03 


0.32 


1.49 


0.03 


13 


E 


0.1335 


3.538 


0.306x102 


0.330 


0.82 


0.02 


0.25 


0.91 


0.04 


14 


E 


0.1664 


4.051 


0.191x102 


0.321 


0.54 


0.04 


0.19 


1.32 


0.03 


15 


E 


0.2107 


4.702 


0.108x102 


0.298 


0.57 


0.05 


0.15 


1.02 


0.05 


16 


E 


0.2726 


5.685 


0.520x10 


0.262 


0.66 


0.04 


0.12 


2.03 


0.00 


17 


E 


0.3679 


9.154 


0.101x10 


0.189 


0.69 


0.05 


0.25 


1.40 


0.02 


9 


F 


0.0603 


2.319 


0.122x10^ 


0.324 


0.73 


2.40 


1.15 


0.77 


0.06 


10 


F 


0.0734 


2.806 


0.705x102 


0.329 


0.65 


2.17 


0.87 


0.48 


0.03 


11 


F 


0.0894 


3.341 


0.421x102 


0.331 


0.65 


1.58 


0.64 


0.66 


0.07 


12 


F 


0.1094 


3.946 


0.255x102 


0.331 


0.58 


0.96 


0.51 


1.03 


0.07 


13 


F 


0.1347 


4.618 


0.153x102 


0.322 


0.55 


0.41 


0.40 


0.56 


0.06 


14 


F 


0.1679 


5.344 


0.930x10 


0.311 


0.56 


0.11 


0.30 


1.45 


0.02 


15 


F 


0.2128 


6.266 


0.524x10 


0.292 


0.63 


0.02 


0.24 


2.05 


0.02 


16 


F 


0.2759 


7.501 


0.251x10 


0.249 


0.57 


0.04 


0.22 


0.85 


0.03 



Table 6: - continued. 
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X bin 


(5 bm 


(x) 


{Q ) 


<i, / 75 

a ja^ 


^stat. 


^rad. 


S-model 








[GeV J 




[%] 


[07' 1 

[%] 


[%] 


1 


A 


0.0069 


0.291 


0.992 


2.19 


1.52 


3.85 


2 


A 


0.0102 


0.404 


0.952 


1.93 


0.94 


2.53 


3 


A 


0.0142 


0.514 


0.979 


1.27 


0.43 


2.01 


4 


A 


0.0186 


0.626 


0.984 


1.41 


0.39 


1.71 


5 


A 


0.0244 


0.725 


0.977 


1.21 


0.10 


1.56 


6 


A 


0.0325 


0.872 


0.965 


1.26 


0.57 


1.47 


7 


A 


0.0398 


0.903 


0.972 


1.94 


0.82 


1.62 


8 


A 


0.0486 


0.929 


0.962 


3.24 


0.54 


2.70 


9 


A 


0.0595 


1.102 


0.948 


3.72 


1.49 


1.37 


10 


A 


0.0724 


1.127 


0.936 


3.69 


0.63 


1.17 


11 


A 


0.0882 


1.167 


0.954 


3.51 


1.84 


0.89 


12 


A 


0.1078 


1.205 


0.923 


5.02 


1.51 


0.81 


13 


A 


0.1330 


1.263 


0.885 


4.14 


1.97 


0.86 


14 


A 


0.1657 


1.337 


0.927 


3.57 


2.48 


1.09 


15 


A 


0.2093 


1.447 


0.873 


2.87 


1.24 


1.90 


16 


A 


0.2642 


1.749 


0.848 


1.72 


1.09 


2.76 


17 


A 


0.3477 


2.672 


0.798 


1.88 


0.09 


2.08 


18 


A 


0.4605 


4.224 


0.783 


1.59 


0.19 


1.77 


19 


A 


0.6346 


8.688 


0.743 


3.90 


0.07 


1.28 


1 


B 


0.0085 


0.354 


0.968 


1.97 


0.69 


2.97 


2 


B 


0.0117 


0.477 


0.967 


1.17 


0.99 


2.49 


3 


B 


0.0152 


0.612 


0.972 


0.93 


0.81 


2.28 


4 


B 


0.0193 


0.770 


0.969 


0.97 


0.33 


2.07 


5 


B 


0.0250 


0.896 


0.981 


0.78 


0.27 


1.69 


6 


B 


0.0327 


1.144 


0.972 


1.00 


0.45 


1.60 


7 


B 


0.0401 


1.074 


0.964 


1.56 


0.32 


1.38 


8 


B 


0.0490 


1.123 


0.962 


1.55 


0.07 


1.30 


9 


B 


0.0596 


1.281 


0.969 


2.96 


0.57 


0.87 


10 


B 


0.0725 


1.359 


0.928 


2.63 


0.18 


0.61 


11 


B 


0.0884 


1.399 


0.925 


2.71 


0.14 


0.55 


12 


B 


0.1082 


1.554 


0.921 


1.73 


0.36 


0.50 


13 


B 


0.1333 


1.632 


0.901 


2.57 


0.73 


0.53 


14 


B 


0.1661 


1.811 


0.903 


1.45 


0.17 


0.54 


15 


B 


0.2104 


2.050 


0.885 


2.13 


0.25 


0.61 


16 


B 


0.2722 


2.471 


0.860 


1.43 


0.06 


0.82 


17 


B 


0.3620 


3.790 


0.809 


1.21 


0.09 


0.60 


18 


B 


0.4900 


6.026 


0.768 


0.91 


0.03 


0.66 


19 


B 


0.6641 


12.778 


0.756 


3.22 


0.09 


0.67 


5 


C 


0.0270 


1.092 


0.990 


1.04 


0.43 


2.03 



Table 7: Results on the inelastic Born cross-section ratio a'^/aP. The statistical uncertainty 
^stat. , the systematic uncertainty Srad. due to radiative corrections and Smodei due to the model 
dependence outside the acceptance are given in percent. Corresponding x bin numbers and 
bin numbers and the average values of x and are listed in the first four columns. The overall 
normalization uncertainty is 1.4%. 
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X bin 


bm 


{x) 


1 r\1\ 

{Q) 




^stat. 


^rad. 


^model 








T 7-21 

[GeV J 




[%\ 


[%] 


[%] 


6 


C 


0.0338 


1.391 


0.990 


1.49 


0.73 


2.15 


7 


C 


0.0401 


1.293 


0.975 


1.18 


0.07 


1.45 


8 


C 


0.0490 


1.460 


0.952 


1.02 


0.19 


1.12 


9 


C 


0.0596 


1.466 


0.941 


1.59 


0.20 


0.75 


10 


C 


0.0726 


1.627 


0.943 


1.59 


0.26 


0.53 


11 


c 


0.0885 


1.691 


0.945 


1.60 


0.08 


0.46 


12 


c 


0.1082 


2.013 


0.917 


1.39 


0.05 


0.39 


13 


c 


0.1334 


2.158 


0.920 


1.57 


0.17 


0.36 


14 


c 


0.1662 


2.458 


0.891 


1.74 


0.10 


0.33 


15 


c 


0.2105 


2.851 


0.862 


1.07 


0.01 


0.31 


16 


c 


0.2724 


3.479 


0.851 


0.96 


0.08 


0.30 


17 


c 


0.3622 


5.125 


0.813 


0.87 


0.03 


0.22 


18 


c 


0.4957 


8.280 


0.760 


0.97 


0.02 


0.23 


7 


D 


0.0408 


1.615 


0.969 


1.01 


0.51 


1.94 


8 


D 


0.0496 


1.927 


0.978 


0.90 


0.24 


1.63 


9 


D 


0.0596 


1.675 


0.961 


1.76 


0.22 


0.74 


10 


D 


0.0726 


1.928 


0.951 


2.11 


0.12 


0.51 


11 


D 


0.0885 


2.080 


0.946 


1.31 


0.10 


0.43 


12 


D 


0.1083 


2.490 


0.914 


1.46 


0.31 


0.34 


13 


D 


0.1334 


2.795 


0.923 


1.20 


0.21 


0.29 


14 


D 


0.1663 


3.195 


0.896 


1.29 


0.10 


0.25 


15 


D 


0.2106 


3.709 


0.877 


0.99 


0.10 


0.21 


16 


D 


0.2726 


4.527 


0.834 


1.25 


0.02 


0.18 


17 


D 


0.3624 


6.797 


0.808 


0.82 


0.02 


0.13 


18 


D 


0.5038 


11.344 


0.770 


1.43 


0.02 


0.31 


9 


E 


0.0597 


1.906 


0.950 


1.31 


0.18 


0.81 


10 


E 


0.0726 


2.276 


0.954 


1.34 


0.20 


0.58 


11 


E 


0.0885 


2.587 


0.944 


1.00 


0.14 


0.44 


12 


E 


0.1083 


3.071 


0.940 


0.81 


0.04 


0.35 


13 


E 


0.1335 


3.538 


0.912 


1.20 


0.04 


0.27 


14 


E 


0.1664 


4.051 


0.901 


0.79 


0.06 


0.21 


15 


E 


0.2107 


4.702 


0.878 


0.82 


0.04 


0.17 


16 


E 


0.2726 


5.685 


0.853 


0.95 


0.04 


0.13 


17 


E 


0.3679 


9.154 


0.802 


0.98 


0.02 


0.26 


9 


F 


0.0603 


2.319 


0.969 


1.07 


0.14 


1.26 


10 


F 


0.0734 


2.806 


0.955 


0.96 


0.15 


0.94 


11 


F 


0.0894 


3.341 


0.936 


0.96 


0.09 


0.72 


12 


F 


0.1094 


3.946 


0.947 


0.85 


0.08 


0.56 


13 


F 


0.1347 


4.618 


0.920 


0.81 


0.05 


0.44 


14 


F 


0.1679 


5.344 


0.898 


0.81 


0.03 


0.32 


15 


F 


0.2128 


6.266 


0.874 


0.91 


0.01 


0.26 


16 


F 


0.2759 


7.501 


0.829 


0.82 


0.04 


0.23 



Table 7: - continued. 
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Parameter 


ATT 1\T07 T3 


GD07-P 


GDll-P 


GDll-D 




value 


value 


unc. 


value 


unc. 


value 


unc. 


™2 7'2i 
mg[LreV J 


u.oiyoo 


0.454 


0.0283 


0.5063 


0.0236 


0.426 


- 


m^lKjev \ 


/in /I r^v 
4y.40 / 


30.7 


2.85 


34.75 


2.56 


0.00007713 


0.000319 




U.ioUoz 


0.117 


0.0465 


0.03190 


0.00828 


0.2293 


0.0763 




U.0Z044 


1.15 


0.358 


1.374 


0.349 


2.65 


- 


A 2 rr'^A^2i 


U.UOOZ 1 


0.06527 


- 


0.06527 


- 


0.06527 


- 


a-pi 


-U.UoUo 


-0.105 


0.00507 


-0.11895 


0.00411 


-0.4287 


0.0773 


a-p2 


n /I /I Q 1 
U.44oiZ 


-0.495 


0.0306 


-0.4783 


0.0224 


0.2891 


0.0646 


a-p3 


i.i ^ uy 


1.29 


0.243 


1.353 


0.203 


0.3931 


0.0935 


0-P4 


U.ODzyz 


-1.42 


0.491 


1.0833 


0.0908 


-27.212 


0.649 


0-P5 


1 oni ^ 
i.oyi r 


4.51 


0.540 


2.656 


0.363 


30.687 


0.650 


0-P6 


1 Q /I on 


0.551 


0.120 


1.771 


0.204 


0.04577 


0.00198 


C-P7 


n ocnf^'y 
U.zoUo / 


0.339 


0.0194 


0.3638 


0.0140 


0.00073 


0.00181 


C-P8 


n oooni 


n 1 97 




n 1 91 1 


n ni fi'^ 


n Q7zl1 


u.uouu 


C'PQ 


2.iy7y 


1.16 


0.246 


1.166 


0.197 


0.8722 


0.0800 




0.584 


0.374 


0.0320 


0.3425 


0.0243 


0.2986 


0.0237 


«-R,2 


0.37888 


0.998 


0.102 


1.0603 


0.0640 


3.615 


0.230 


«-R,3 


2.6063 


0.775 


0.112 


0.5164 


0.0440 


1.1455 


0.0653 


bm 


0.01147 


2.71 


0.393 


-10.408 


0.627 


1.987 


0.333 


bizb 


3.7582 


1.83 


0.537 


14.857 


0.627 


7.150 


0.343 


bne 


0.4y338 


1.26 


0.296 


0.07739 


0.00861 


0.9350 


0.0935 


cm 


0.80107 


0.838 


0.106 


1.3633 


0.0867 


1.0316 


0.0783 


ens 


0.97307 


2.36 


0.557 


2.256 


0.503 


26.36 


6.19 


C7^9 


3.4942 


1.77 


0.209 


2.209 


0.260 


3.024 


0.136 



Table 8: Results of fits to the cross section (Jl+t using the ALLM functional form Parameters 
of the ALLM97 fit and the GD07-P fit (both based on proton data) arc quoted together with 
the new parameters determined from current world data including the HERMES results discussed in 
this paper on proton (GDll-P) and deuteron (GDll-D) cross sections. In Ref. the uncertainties 
were calculated with an UP value equal to 24.7. That identical fit was redone here with the same 
value of UP=1 as the new fits in order to provide an appropriate comparison of its uncertainties. 
The ALLM97 fit provided no uncertainties. The uncertainties given in the table for GD07-P, GDll- 
P, and GDll-D correspond to the diagonal elements of the full covariance matrix which must be 
used to calculate uncertainties in F2 or cross sections. They can be obtained by contacting the 
management of the HERMES collaboration: <management@hermes. desy.de>. The parameter Ag 
has no uncertainty as it is fixed in the fits. In the case of the deuteron fit, the parameters rrig and 
Qq had to be fixed in order to obtain a good fit. 
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Data set 



GDll-P 



GDll-D 



HERA (positron beam) 
HERA (electron beam) 
E665 m 
NMC-90 GeV 
NMC-120 GeV 
NMC-200 GeV 111 



NMC-280 GeV H 



BCDMS-100 GeV 
BCDMS-120 GeV 
BCDMS-200 GeV 
BCDMS-280 GeV 



SLAG E49a [|T9| 
SLAG E49b g9| 
SLAG E61 



19 



SLAG E87 m 
SLAG E89a pl 



SLAG E89b |19| 



SLAG E139 |19| 



SLAG E140 PI 
JLAB 



JLAB [2C] 



JLAB |21|] 

JLAB (Rosenblutli) 12 
JLAB (Model Dependent) U 
HERMES (T.A.) 



-0.0065 
-0.0067 
0.020 
-0.00020 
0.011 
0.0093 
0.0035 
-0.032 
-0.028 
-0.027 
-0.023 
0.016 
0.022 
0.016 
0.016 
0.036 
0.018 



-0.012 
-0.0063 

0.0014 
0.0085 
0.015 



-0.014 
-0.029 
-0.0096 
0.0036 
0.0023 

-0.0075 
-0.0060 
-0.0031 
-0.0013 
0.0062 
0.0070 
0.0045 
0.0087 
0.00081 
0.0014 
0.0025 
-0.040 

-0.0012 
0.0088 
0.0088 
-0.022 



Table 9: Normalizations Ski^k (sec Eq. (6.2)) obtained in the GDll fits. 
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